Non-local transport in normal-metal/superconductor hybrid structures: the role of 

interference and interaction 
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We have measured local and non-local conductance of mesoscopic normal-metal/superconductor 
hybrid structures fabricated by e-beam lithography and shadow evaporation. The sample geometry 
consists of a superconducting aluminum bar with two normal-metal wires forming tunnel contacts 
to the aluminum at distances of the order of the superconducting coherence length. We observe 
subgap anomalies in both local and non-local conductance that quickly decay with magnetic field 
and temperature. For the non-local conductance both positive and negative signs are found as a 
function of bias conditions, indicating at a competition of crossed Andreev reflection and elastic 
cotunneling. Our data suggest that the signals are caused by a phase-coherent enhancement of 
transport rather than dynamical Coulomb blockade. 
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I. INTRODUCTION 



Andreev reflection^ (AR) is the process responsible for 
the transfer of electrons from a normal metal into a su- 
perconductor at energies below the superconducting gap 
energy. In multi-terminal structures with two (or more) 
normal metals (N) connected to a single superconduc- 
tor (S), non-local or crossed Andreev reflection (CAR) 
may occur, where an electron entering the superconduc- 
tor from one normal-metal contact A forms a Cooper 
pair by emitting a hole into a second contact A 
competing process is elastic cotunneling (EC),'* where 
an electron is transmitted to contact B without the for- 
mation of a Cooper pair. CAR has attracted theoretical 
and experimental^i^i^i^iiSiiiiii^iii^ attention mainly be- 
cause it is predicted to create spatially separated, entan- 
gled electrons in a solid-state environment (see Ref. [13 
for a brief review). While CAR can be readily probed 
by spin selection^i^ using ferromagnetic electrodes, this 
approach is unsuitable for entangler devices, since pro- 
jecting the spin will destroy entanglement. Therefore, 
understanding the competition between CAR and EC in 
multiterminal NSN structures is a prerequisite for the 
successful design of superconducting solid-state entan- 
glers. In a recent experiment a bias-dependent crossover 
from EC to CAR has been observed in a diffusive 
NSN structure with low-transparency tunnel contactsi^ 
While both quantum mechanical interference^ii^ii^iiiii^ 
and Coulomb interactio n*^!*^ have been proposed to ex- 
plain the result, its origin is not yet clear. Here, we 
present an experimental investigation of both local and 
non-local transport in lateral NSN hybrid structures. 
We observe a subgap anomaly in the non-local conduc- 
tance similar to the observations in Ref. 0- Comparison 
with local-transport data reveals that the non-local sig- 
nal is controlled by a phase-coherent enhancement of lo- 
cal transport. Dynamical Coulomb blockade is shown to 
be present in the samples, but can be clearly ruled out 



as the cause of the non-local conductance signal by its 
different dependence on temperature and magnetic field. 



II. EXPERIMENTAL 

Figure [T]^a) shows an SEM image of the layout of our 
samples. They mainly consist of a narrow superconduct- 
ing aluminum wire of about 60 /im length (contacts 1, 2 
and 8), and perpendicular to it, two copper wires (con- 
tacts 3/4 and 5/6) forming tunnel junctions A and B to 
the aluminum. A third tunnel junction (contact 7) was 
used for some control experiments. The central junction 
area is shown on an enlarged scale in panel (b), together 
with a typical measurement configuration. Three metal 
layers were evaporated onto an oxidized silicon substrate 
under different angles through a shadow mask fabricated 
by standard e-beam lithography technique. First, an 
auxiliary layer of 30 nm thick copper (Cul), was evap- 
orated to form Ohmic contacts to the two subsequent 
layers at interconnects to the lines 2,3,6,7 and 8. In the 
second evaporation step under a different angle, the su- 
perconductor strip, an aluminum bar of 30 nm thickness, 
was created. The aluminum film was then oxidized in 
situ by applying « 1 Pa of pure oxygen for « 5 min to 
form a tunnel barrier. Then the third layer (Cu2), made 
of copper with thickness tcu2, was evaporated, forming 
the two tunnel contacts A and B to the aluminum. Sam- 
ple parameters such as contact areas, resistivities, etc., 
varied slightly between fabrication batches. Consistent 
results were obtained from six different samples, with 
parameters given in Table H] We show here mainly data 
from samples I and II. 

The samples were mounted into a shielded box ther- 
mally anchored to the mixing chamber of a dilution re- 
frigerator. The measurement lines were fed through a 
series of filters to eliminate RF and microwave radiation 
from the shielded box. A voltage Vex consisting of a vari- 
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FIG. 1: (Color online) (a) SEM image of the sample layout, 
(b) Close-up of the contact region, colorized for clarity. Two 
copper (orange) wires (A and B) are connected by overlap 
tunnel junctions to a weakly oxidized aluminum bar (blue). 
The copper wires fork near the tunnel contact to allow four- 
probe characterization. As a by-product of shadow evapora- 
tion, additional unconnected parts of the sample (grey), and 
a Cu/Al/oxide/Cu trilayer (green), are formed. The contact 
configuration and measurement scheme is also indicated. The 
three metal layers (Cu2, Cul, Al) as created by shadow evap- 
oration are indicated above the image. 



TABLE I: Characteristic parameters of the six samples. Film 
thickness tcu2, contact distance d, coherence length ^, alu- 
minum resistivity pAi at T = 4.2 K, normal-state tunnel junc- 
tion conductances Ga and Gb, and contact transparencies 7a 
and Tb- 



sample 


tCu2 


d 




PAl 


Ga Gb 


Ta Tb 




(nm) 


(nm) 


(nm) 


(/iOcm) 


(mS) 


(xl0~5) 


I 


20 


75 


120 


4.1 


1105 1530 


4.9 5.6 


II 


20 


100 


120 


4.2 


660 760 


5.3 5.2 


III 


20 


115 


170 


2.1 


520 380 


4.6 3.1 


IV 


30 


300 


140 


3.1 


1020 1150 


3.3 3.3 


V 


30 


105 


110 


5.1 


1480 1520 


7.7 7.4 


VI 


30 


105 


110 


5.3 


1390 1110 


7.3 5.6 



able DC bias and a low-frequency AC excitation was ap- 
plied to contact A, and the resulting currents /a and 
/b through both contacts were measured using indepen- 
dent current amplifiers. The actual voltage Va across 
the contact was measured via additional leads in a four- 
probe configuration. Interchanging current and voltage 
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FIG. 2: (Color online) X-ray diffraction patterns of an un- 
patterend Al/oxide/Cu sandwich, with substrate background 
subtracted, (a) 6 — 26 scan, with expected positions and rela- 
tive intensities of (111) and (200) reflections for both Al and 
Cu indicated by vertical bars, (b) Rocking curve with angle 
of inclination u) for the (111) reflections of Al and Cu. 



probes did not change the observed signals. Measure- 
ments were also performed with the roles of injector and 
detector interchanged between contact A and B, with 
consistent results. In some cases, an additional DC bias 
V-Q was applied to contact B. Voltage and current po- 
larities are indicated in Fig. [T]^b) by plus signs and ar- 
rows, respectively. For our choice of current polarities 
EC leads to a positive non-local conductance, while CAR 
leads to a negative signal. The AC components of Va, 
/a and /b were measured simultaneously with lock-in 
technique, and from the in-phase signals the local and 
non-local differential conductances 5aa — dI\/dVf^ and 
5AB = diB/dVA were extracted. Since we are interested 
in non-local conductances which are smaller than local 
conductances by orders of magnitude, care was taken 
to ensure that the measured signals were not affected 
by phase shifts or crosstalk between measurement lines. 
We performed extensive electronic circuit simulations of 
the entire measurement setup including amplifiers, fil- 
ters, and cryostat wiring. The simulations showed that 
the signals are reliable for measurement frequencies up 
to / « 200 Hz, which was also confirmed experimentally. 
All data shown in this paper were taken at / w 37 Hz 
with a typical AC amplitude of 5 /iV. In the remainder 
of the paper, Va and Vb will refer to DC bias voltage, G 
will refer to normal-state conductance, and g will refer 
to differential conductance. 

For X-ray characterization, we also prepared unpat- 
terned bilayers in the same way as the overlap contact 
area of the patterned samples, i.e., with a first layer of 
aluminum, subsequent oxidation, and a second copper 
layer on top. X-ray diffraction of these bilayer films was 
performed using copper Ka radiation and a solid-state 
detector. Figure [2Ja) shows 9 — 26 scans, with substrate 
background subtracted. The only reflections we observed 
were the (111) reflections for both Al and Cu. As an ex- 
ample, the expected positions of the (200) reflections, and 
their expected intensities relative to the (111) reflections, 
are also indicated by vertical bars for both materials. In 
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FIG. 3: (Color online) Local conductance qaa (a,c) and non- 
local conductance qab (b,d) of sample I as a function of tem- 
perature T at Va = (a,b) and injector bias Va at T = 25 mK 
(c,d). The insets in (a) and (b) show the low-temperature re- 
gion of both plots on an enlarged scale. The solid lines in (a) 
and (c) are fits to a standard BCS tunneling model. 



FIG. 4: (Color online) Local differential conductance qaa 
(a,c) and non-local differential conductance qab (b,d) as a 
function of injector bias Va taken at different temperatures T 
for sample I (a,b), and different magnetic fields B for sample 
II (c,d). 



Fig.m^b), rocking curves for the two (111) reflections are 
shown, with an FWHM of 0.35° and 0.26° for Al and 
Cu, respectively. From the absence of all but the (111) 
reflections, and their small width in the rocking curve, 
we conclude that the films have a nearly perfect (111) 
texture. 



III. RESULTS 

The local conductance ^aa of the injector contact A 
of sample I is shown in Fig. [3Ia) as a function of tem- 
perature T for Va = 0. Above the critical temperature 
Tc = 1.35 K of the aluminum film, a constant normal- 
state tunnel conductance is observed. Below Tc, the 
conductance gradually drops to zero, and can be fitted 
with a standard BCS^-'^^ tunneling model, with zero- 
temperature gap Aq = 205 /ieV and normal-state con- 
ductance Ga = 1105 /iS. Below about 250 mK, an addi- 
tional conductance contribution not covered by the tun- 
neling model is observed. This anomaly can be better 
resolved in the inset, where it is shown on an enlarged 
scale. Figure ^h) shows the non-local conductance ^ab 
measured simultaneously. The overall behavior is sim- 
ilar to the local conductance, including the presence of 
a low-temperature anomaly, as seen in the inset. The 
low-temperature anomaly decreases more steeply with in- 
creasing temperature as compared to the local anomaly 
seen in the inset of Fig. [3Ua), as will be discussed later. 



The local differential conductance gAA for the same 
configuration is shown in Fig. [3ljc) as a function of bias 
Va at T = 20 mK. It has the form of the BCS den- 
sity of states as expected for a tunnel contact with low 
transparency (T « 6 • 10^^ for this sample, estimated 
from Ga and the contact area). The solid line is a fit 
to the BCS tunneling model, with the same parameters 
as in Fig. Efa). Except for the small zero-bias, low- 
temperature anomaly already seen in panel (a), and dis- 
cussed further below, the subgap conductance is negligi- 
ble, indicating a high quality tunnel barrier without pin 
holes. The non-local differential conductance gAB corre- 
sponding to panel (c) is shown in panel (d). Similar to 
the local signal, the subgap conductance is almost zero, 
with a sharp transition to a flnite signal above the gap. 
The signal above the gap can be attributed to charge im- 
balance caused by quasiparticle injection. A detailed in- 
vestigation of the non-local charge-imbalance signal will 
be reported elsewhere. 

Figure [5] shows the low-temperature subgap anoma- 
lies seen in both local and non-local conductance on an 
enlarged scale. Panel (a) shows the local differential con- 
ductance of sample I for different temperatures in the 
range between 25 mK and 200 mK, at zero magnetic 
fleld. A sharp peak is observed at zero bias. The peak 
height decreases with increasing temperature, as already 
observed in the inset of Fig.ISja). At Va ~ ±20 /.tV, min- 
ima are seen, and for the lowest temperatures the differ- 
ential conductance actually becomes negative. At higher 
bias, there are a series of side-maxima, which will be dis- 
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Vb (uv) (uV) 
FIG. 5: (Color online) 

(a) gAB as a function of Va for different Vb ■ 

(b) Oscillatory part Qosc of the non-local conductance qab as 
a function of both Va and Vb- 

(c) Qosc at Va = (symbols) and local conductance gBB of 
contact B (line) as a function of Vb- 

(d) pshift (symbols) and derivative of the local conductance 
dgBB /rfVe as a function of Vb (line) . 

All data are taken at B = and T — 20 mK from sample II. 
For an explanation of Qosc and gshUt, see text. 

cussed in more detail below (see Fig. |8]). The positions 
of the maxima and minima are independent of tempera- 
ture. Panel (b) shows the non-local differential conduc- 
tance measured simultaneously. The non-local conduc- 
tance shows a structure similar to the local one, with a 
central peak and negative side-minima. The minima oc- 
cur at Va ~ ±15 fiV, sligthly below the position of the 
minima in the local conductance. The anomaly decreases 
with temperature, and drops below the noise floor above 
100 niK. Panels (c) and (d) show the local and non-local 
differential conductances of sample II at low tempera- 
ture for different magnetic fields B applied parallel to the 
substrate plane along the direction of the copper wires 
forming the contacts. The signals are almost identical to 
those observed in sample I. As temperature, the magnetic 
field leads to a suppression of the anomalies, on a field 
scale much below the critical field of the aluminum wire 
(Be ~ 600 mT, determined from conductance measure- 
ments) . While the amplitude of the signals is suppressed 
by both finite temperature and magnetic field, no change 
in the overall shape or voltage scales is seen. 

To study the non-local conductance pattern in more 
detail, we have applied an additional DC bias Vb to the 
second contact. Figure El^a) shows ^ab as a function of 
Va for three different values of Vb- The effect of Vb is 
twofold: first, there is an overall vertical shift indepen- 



FIG. 6: (Color online) Normalized amplitudes of gosc and 
ffshift, together with the square of the amplitude of the subgap 
anomaly in the local conductance (5aa)^, as a function of 
temperature T (a) and applied magnetic field B (b). 
Data taken from sample I. 

dent of Va, as can be clearly seen at large bias. Second, 
there is a horizontal shift of the maxima and minima of 
the oscillatory pattern. We can therefore describe the 
signal as 

ffAB(VA, Vb) = gshift(VB) + 5osc(Va, Vb) 

The vertical shift ^shift can be extracted from the data 
at bias voltages around 60 /iV, where the oscillatory sig- 
nal (7osc has died out. gosc obtained by subtracting ^shift 
is shown in Fig. [5jb) in a color plot as a function of both 
Va and Vb . The pattern is limited to a bias voltage win- 
dow |Va, Vb| < 40 /iV, and the shift of the maxima and 
minima is linear in Vb, with a slope of about 0.7 (i.e., the 
signal does not simply depend on Va — Vb). It can also 
be seen that additional maxima and minima shift into 
the bias voltage window upon increasing |Vb|. 

Figure [SJc) shows gosc as a function of Vb for Va = 
(i.e., a vertical cut through Fig. [5jb), together with the 
local conductance ^bb of contact B, measured indepen- 
dently. As can be seen, the two signals roughly scale onto 
each other. A further correlation between non-local and 
local data can be seen in Fig. E^d), where ^shift is plot- 
ted as a function of Vb , along with the second derivative 
dgBB/dVB — (PlB/dV^ of the local conductance of con- 
tact B. 

Figure [H^a) and (b) show the dependence of the am- 
plitudes of the local and non-local subgap anomalies as 
a function of temperature and magnetic field. Data 
are taken from sample I, and normalized to the low- 
temperature and low-field values. Both the local and 
non-local signals decay rapidly, at temperature and field 
scales well below the critical temperature Tc — 1.3 K 
and critical field Be — 600 mT, respectively. As can be 
seen, the two components gosc and gshift of the non-local 
conductance both scale with the square of the amplitude 
of the zero-bias anomaly in the local conductance. The 
steeper decrease of the non-local conductance as a func- 
tion of temperature has already been noted in Fig. [31 

In order to elucidate the role of dynamical Coulomb 
blockade (DCB) in our samples, we have also investigated 
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FIG. 7: (Color online) Local conductance qaa (a) and non- 
local conductance gAB (b) of sample II in the normal state 
at B = 1.5 T for different temperatures T. The solid line in 
panel (a) is a fit to the standard model of dynamical Coulomb 
blockade.— 



the differential conductance in the normal state, where 
the impact of DCB is well known. Figure [7] shows both 
local (a) and non-local (b) differential conductance in 
a magnetic field of 1.5 T, sufficiently large to suppress 
superconductivity, at different temperatures. At lowest 
temperature, the local conductance shows a low-bias dip 
which can be well described by DCB in the presence of 
an Ohmic environment,^^ with environmental impedance 
R^^^ = 160 see fit in Fig. [T^a). The Coulomb dip 
persists up to temperatures of about 2 K, well above the 
critical temperature of aluminum. A similar dip appears 
in the non-local conductance, shown in panel (b). The 
dip in the non-local conductance is slightly narrower than 
in the local conductance, and also persists up to about 
2 K. 



IV. DISCUSSION 

A. Local conductance 

We will first focus on the interpretation of the subgap 
anomaly of the local conductance. Figure [Hl^a) shows the 
differential conductance ^aa of sample II at low temper- 
ature and zero field, together with three different models 
described in the following. Here, a small regular back- 
ground contribution of about 1 /iS has been subtracted 
to show the anomaly more clearly. Panel (b) shows the 
same data on an enlarged scale, together with the ballis- 
tic model. 

Diffusive model: A zero-bias conductance peak is 
known to occur in low-transparency tunnel junctions 
between diffusive normal metals and superconductors 
as a result of phase-coherent enhancement of Andreev 
reflection (reflectionless tunneling).— i^liS^ The subgap 
conductance due to reflectionless tunneling for differ- 
ent experimental situations has been calculated in Refs. 
[26,27,28,29. A fit to the model of an overlap junction 
including pair-breaking (equation (6) of Ref. ^9) shows 
good agreement for the central peak, but fails to describe 



FIG. 8: (Color online) (a,b) Differential subgap conduc- 
tance of sample II at _B = and T = 20 mK, with reg- 
ular background subtracted (symbols). The lines are three 
different fits described in detail in the text, (c) Measured 
non-local conductance amplitude qab vs. model prediction 
<7ab = ■ Qaa ■ Qbb ■ exp(— d/^). See text for details. 



the addtional minima and maxima. It should be noted 
that the characteristic energy scale en = 0.18 /ieV con- 
trolling the weight of the anomaly was actually not fitted 
but simply calculated from known sample parameters. 
Only a small amount of pair-breaking (7n — 4 ^eV) had 
to be adjusted to obtain the correct width (see Ref. [2^ 
for the definition of the parameters). 

Ballistic model: Reflectionless tunneling results from 
the constructive interference of repeated attempts of an 
electron to be Andreev-reflected at the tunnel junction. 
In a thin-film overlap junction these repeated attempts 
will mostly come from electrons bouncing back and forth 
between the junction and the upper surface of the normal 
metal film. Elastic scattering in thin films can largely be 
attributed to grain boundaries and surfaces, since both 
the elastic mean free path and the grain size are usually 
of the order of the film thickness. Therefore, while elec- 
tron motion in the lateral direction is certainly diffusive, 
it can be assumed to be close to ballistic perpendicular to 
the films, especially since our films show a nearly perfect 
(111) texture. The ballistic limit of reflectionless tun- 
neling has been considered in Ref. [10 using semiclassical 
trajectories, similar to the treatment of diffusive motion 
in Ref. [23- Including the phase difference due to flnite 
electron energy E, Acj) = 2EL/'hvp^ where L is the 
length of a trajectory and up is the Fermi velocity, we ob- 
tain a differential conductance spectrum with a central 
peak and small, equidistant side-maxima, as seen best 
on a larger scale in Fig. [5]^b). The flt parameters were 
chosen such that the height of the central peak and the 
position of the side-maxima are reproduced. Better fits 
for either the central peak or the side-maxima could be 
obtained, but only at the expense of the other. In addi- 
tion to the position of the side-maxima, also their decay 
at higher bias is reproduced. The negative differential 
conductance remains unexplained in this model. 

Josephson model: A zero-bias peak with a subsequent 
region of negative differential conductance can be ob- 
served in voltage biased Josephson junctions in the phase- 



6 



diffusion regime. A part of the normal-metal leads of our 
junctions, the green shaded region in Fig. [TJb), actu- 
ally consists of a Cu/Al/oxide/Cu trilayer as a result of 
shadow evaporation. While superconductivity in the alu- 
minum of this trilayer must be strongly suppressed by the 
inverse proximity effect, we can not exclude the presence 
of a weak Josephson coupling between the trilayer and 
the aluminum bar. As in the case of reflectionless tunnel- 
ing, the bottleneck for the Josephson coupling would be 
pair transmission through the tunnel barrier. The model 
of phase diffusion^^ plotted in Fig. [5Ja) gives a quali- 
tative description of the central peak and negative dif- 
ferential conductance, but no side-maxima are predicted 
here. The phase diffusion model predicts that the current 
maximum (i.e., the zero of the differential conductance) 
appears at a voltage which is proportional to temper- 
ature. In contrast, no temperature dependence of the 
voltage scales is observed in our samples. In the light 
of the Josephson model, one may ask whether the ad- 
ditional side-maxima are the result of multiple Andreev 
reflection (MAR)i'^^ Peaks due to MAR would appear at 
characteristic voltages V = 2A/ne, where n is an integer 
corresponding to the correlated transfer of n electrons. 
MARs are inconsistent with our data for several reasons: 
first, they are observable only at relatively high contact 
transparency due to the multiple transmission through 
the interface. For the same reason, the features with 
smallest n, i.e., at highest bias, should be the largest, 
in contrast to our observations. Also, the observed side- 
maxima are equidistant, and too few for MAR (six should 
be visible between 50 and 150 /ieV, assuming the bulk 
gap A = 200 (leY). 

All three models have in common that they require 
phase-coherent motion of Andreev pairs on the normal 
metal side of the tunnel junction. This is consistent with 
the fast suppression of the conductance anomaly both 
as a function of temperature, and in the presence of a 
magnetic field. We conclude that the oscillatory struc- 
ture in the local conductance is caused by reflectionless 
tunneling close to the ballistic limit, while the negative 
differential conductance might indicate the presence of a 
weak Josephson coupling. 

Finally, we would like to address the impact of dynam- 
ical Coulomb blockade (DCB) in the presence of a finite- 
impedance electromagnetic environment on the local con- 
ductance. As already discussed (see Fig. [7]), in the normal 
state the local conductance exhibits a suppression at low 
bias, which can be explained by DCB with an Ohmic 
environment Considering our sample design, we ex- 
pect the series impedance to be mainly given by the long, 
narrow aluminum wire. Indeed, the resistance of the alu- 
minum wire between tunnel junction B and lead 1 (see 
Fig- ID) is 190 il, similar to i?onv — 160 H. extracted from 
the fit. In the superconducting state. Coulomb blockade 
of subgap transport is expected to be stronger than in the 
normal state due to the double charge transfer of Andreev 
refiection. Despite the presence of dynamical Coulomb 
blockade, we observe a phase-coherent enhancement of 



subgap conductance in the superconducting state. 



B. Non-local conductance 

Theoretical descriptions of CAR and EC in lowest or- 
der of contact transparency, i.e., in the tunneling limit, 
predict the probability of both processes to be propor- 
tional to GaGbCj where Ga and Gb are the normal-state 
conductances of the two contacts, and C is a factor related 
to the propagation of virtual quasiparticles between the 
contacts in the superconductor. ( depends on sample 
geometry, the mean free path and the coherence length 
In a one-dimensional geometry with two point con- 
tacts at a distance d, C = exp(— d/^). Since CAR and EC 
contribute with opposite signs to the non-local conduc- 
tance, their effect is expected to cancel in the tunneling 
limit, i.e., g^B = 0. This is in contrast to our observa- 
tion of a positive non-local conductance (corresponding 
to EC) at low bias, followed by a cross-over to a negative 
signal (corresponding to CAR) at higher bias. As can be 
seen in Fig.lHIb), the crossover is actually part of a regu- 
lar oscillatory pattern as a function of both injector and 
detector bias voltage. We have found several correlations 
between the local and non-local subgap anomalies: Both 
exhibit a central peak and an oscillatory dependence on 
bias voltage. Some aspects of the local and non-local sig- 
nals can be scaled onto each other (Figs. (5}; and d). In 
particular, the amplitude of the non-local conductance 
anomaly scales with the square of the local one, both as 
a function of temperature and magnetic field. This shows 
conclusively that both must have a common origin. 

It has been predicted recently that reflectionless tun- 
neling in diffusive NSN structures leads to a subgap 
anomaly in the non-local conductance which is propor- 
tional to the product of the subgap anomalies in the lo- 
cal conductance of injector and detector contact, i.e., 
f/AB oc 5aa5bb'^ This is similar to the prediction of 
the lowest-order tunneling model discussed above, ex- 
cept that the normal-state conductances Ga and Gb are 
replaced by the actual subgap anomalies g^A and .^bb- 
Since the model in Ref. 3^ assumes purely diffusive mo- 
tion, it only predicts a zero-bias peak, but no additional 
oscillatory structure. For one-dimensional systems, the 
quantitative prediction is (^ab ~ -Rj •ffAA-ffBB'exp(— d/^), 
where is the normal-state resistance of the supercon- 
ducting wire on a length of We cannot reproduce the 
non-local conductance pattern in detail with this model, 
and will therefore restrict ourselves to a discussion of the 
signal amplitudes. First we note that the scaling of gAB 
with seen in Fig. [6] is fully consistent with the model, 
since within one sample ^aa ~ ffBB- For a comparison 
of the different samples, we have plotted the measured 
non-local conductance amplitude against the prediction 
in Fig.ISfc). As can be seen, the scaling is obeyed reason- 
ably well, but the observed non-local signals are too large 
by three orders of magnitude. The scaling is expected to 
hold for arbitrary sample geometries, while the quantita- 
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tive prediction is for the specific case of one-dimensional 
systems. In the one-dimensional case, the same model 
predicts the local anomaly to be jaa ~ RnG\, where i?N 
is the resistance of the copper wire over the normal- metal 
coherence length — y^fiD/e. For our sample parame- 
ters, this underestimates the local anomaly by about one 
order of magnitude, whereas the model of an extended 
overlap junction used in Fig. [Sfa) provides the correct 
signal amplitude. Whether a theoretical model of our 
specific non-local geometry would similarly remedy the 
quantitative disagreement seen in Fig. E^c) remains an 
open question. Also, a model treating the impact of re- 
flectionless tunneling on non-local transport in ballistic 
structures would be highly desirable to see whether the 
oscillatory structure seen in the non-local conductance 
can be reproduced. 

An alternative explanation for the finite non-local sub- 
gap conductance is dynamical Coulomb blockade. In 
Ref. [l^, a non-local version of DCB has been predicted 
to lift the exact cancellation of CAR and EC in the 
tunneling limit. Depending on whether the blockade of 
CAR or EC is stronger, the other process dominates, and 
gAB may be either positive or negative. In this model, 
the blockade of CAR and EC is controlled by coupling 
to electromagnetic modes of different symmetry. Since 
these modes, in general, have different energies, a bias- 
dependent crossover from dominating EC to CAR can be 
explained. The situation is simpler in the normal state, 
where only EC is possible. In this case, the expectation is 
a suppression of both local and non-local transport at low 
bias, as indeed observed in Fig.[71 However, the Coulomb 
dip in both local and non-local conductance persists to 



temperatures T > 1 K, which is inconsistent with the 
temperature range T < 150 mK of the subgap anomalies 
in the superconducting state. Also, DCB is not affected 
by a magnetic field (the data in Fig. [7] were taken at 
B = 1.5 T), whereas the subgap anomalies are restricted 
to B < 100 mT. The differences in energy, temperature, 
and magnetic field scales clearly show that Coulomb in- 
teraction, while present in our samples, is not the cause 
of the subgap anomalies in the superconducting state. 



V. CONCLUSION 

In conclusion we have presented experimental re- 
sults on non-local transport in mesoscopic normal- 
metal/superconductor hybrid structures. Local and non- 
local conductance exhibit correlated subgap anomalies, 
which are shown to be due to phase-coherent enhance- 
ment of transport at low energies. A systematic depen- 
dence of non-local transport on bias conditions is ob- 
served, which allows control over the dominating trans- 
port processes. Dynamical Coulomb blockade can be 
ruled out as the cause of the subgap anomalies. A com- 
prehensive theoretical description is highly desirable as 
a guideline for the design of superconducting solid-state 
entanglers. 
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